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Abstract: Since 2002, the number of detected blazars at > 100 GeV has more than doubled. I study 
all currently known BL Lac-type objects with published energy spectra. Their intrinsic energy spectra are 
reconstructed by removing extragalactic background light attenuation effects. The emission properties 
are then compared and correlated among each other, with X-ray data, and with the individual black hole 
masses. In addition, I consider temporal properties of the very high energy 7-ray flux. Key findings 
concern the flux-black hole mass and variability scale-black hole mass connections and the correlation 
of the spectral slope and the luminosity. As a specific application, the study allows to constrain the still 
undetermined redshift of the blazar PG 1553+1 13. 



Introduction and Approach 

In order to assess both the acceleration mecha- 
nisms in blazars and extragalactic background light 
(EBL) absorption effects, not only individual ob- 
ject studies, but also the investigation of a large 
sample of very high energy (VHE) 7-ray emitting 
blazars is desirable. Ideally it should encompass 
a wide range in redshift for EBL studies and also 
include groups of sources at comparable distances 
in order to study intrinsic properties of the indi- 
vidual sources without possible systematic uncer- 
tainties caused by the EBL de-absorption. The 
preconditions for such comparative blazar studies 
have much improved recently. To date, the VHE 
7-ray blazar sample comprises 18 BL Lac objects 
(e.g. [1]), with z = 0.030. . .0.212. 1 A sample 
of all blazars with reported energy spectra is stud- 
ied by inferring the intrinsic emission properties of 
the individual objects and by probing correlations 
of their VHE 7-ray and X-ray emission properties 
with black hole (BH) mass (M,) estimates. 

Black Hole Masses. We estimate M, by evalu- 
ating the M.-cr relation [2], i.e. the tight correla- 
tion of the stellar velocity dispersion a and M, of 
nearby galaxies. This approach assumes that AGN 
host galaxies are similar to inactive galaxies as far 



as the M,-a relation is concerned. We find that the 
currently VHE 7-ray emitting BL Lacs are flatly 
distributed in M. = (10 8 - 10 9 5 )M Q . Note that 
although AGNs harbor BHs with M, > 1O 6 M , 
up to now only BL Lacs with M. > 10 8 M Q have 
been discovered in VHE 7-rays. 

Intrinsic VHE 7-Ray Emission Parameters. 

The photon spectra measured in the VHE range 
suffer absorption by EBL [3], Here, the intrin- 
sic blazar spectra are reconstructed using the EBL 
"low" model given in [4]. For sources that have 
been found in different flux states, "low state" and 
"flare" spectra are considered. Data from Mkn 
421 [5], Mkn 501 [6], 1ES2344 [7], 1ES 1959 
[8], PKS 2155 [9], 1H 1426 [10], PKS 2005 
[11], 1ES1218[12], 1ES2356 & 1ES1101 [13], 
PG1553 [14], Mkn 180 [15], PKS 0548 [16], BL 
Lac [17], and 1ES 1011 [18] have been included. 
Throughout this study, the unknown-redshift ob- 
ject PG 1553+113 is assumed at several possible 
z values, but not included further unless explic- 
itly stated otherwise. The extracted observables are 
the intrinsic luminosity at 500 GeV and the intrin- 
sic photon index T in the region around 500 GeV. 
For both, no extrapolations beyond the spectral fits 



1 . See http://www.mppmu.mpg.de/~rwagner/sources/ 
for an up-to-date list. 
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Figure 1: v 1 L 1 vs. vxLx- The Figure 2: Photon index vs. v 1 L 1 . Additional flare states of 



symbols are explained in Fig. 2. 
The PG 1553 points are for as- 
sumed z — 0.1 and z = 0.3. 



sources are marked by gray circles. The PG 1553 points were 
not included in the fit. Right: Four blazars with low and high 
states. 



are required. The resulting intrinsic photon indices 
vary from Y = 1.5 — 3.3. In most current acceler- 
ation models only Y > 1.5 is allowed. This goes 
in line with indications that the EBL absorption ef- 
fects are still smaller than modeled [13]. 

Correlation of X-ray and 7-Ray Lumi- 
nosity 

In SSC models, the X-ray and the VHE emission 
are closely connected, owing to their common ori- 
gin. Fig. 1 shows v 1 L 1 versus the X-ray luminos- 
ity at 1 keV {vxLx', from Costamante & Ghisellini 
2002). A trend towards a correlation is visible, 
even though a strong correlation might not even be 
expected due to different magnetic fields in the in- 
dividual objects. Note that high thermal contribu- 
tions at 1 keV are unlikely and would imply a very 
high amount of gas and pressure. 

Correlation between Photon Index and 
7-ray Luminosity 

Fig. 2 relates the intrinsic photon indices Y to 
v 1 L 1 . A correlation on the 3.3<r level is found, 
which, within SSC models, is compatible with a 
moving IC peak towards higher energies and an 
IC peak energy < 500 Ge V. Sources with observed 
spectra at individual distinct flux states support this 
correlation. Mkn501 and 1ES2344 show a sim- 
ilar change in spectral slope and a luminosity in- 
crease of A(i/ 7 L 7 ) w 20. The luminosity increase 
of Mkn 421 is much lower with A(^ 7 i 7 ) ss 10. 



Correlations of 7-Ray Emission with 
the BH Properties and z 

The properties of blazar 7-ray emission are ex- 
pected to be connected to BH properties, like M, 
and its spin, since scaling laws govern BH physics, 
in particular length and time scales [20]. Currently, 
only A/, can be reliably estimated. The BH spin 
remains inaccessible by large; the accretion rate 
might be indirectly accessible through the (radio) 
jet power. A first study of the connection of source 
properties and M, of the then-established five TeV 
blazars [21] did not find any correlations, except 
for an indication of a connection between the X- 
ray flare duty cycle and M. (see below). 

Fig. 3 shows the correlation of Y and v 1 L 1 with 
M. and also tests for possible correlations with 
z. The latter are not expected from physics, but 
may identify selection effects in the data sample 
and/or an inaccurate EBL model. Only sources 
with hard intrinsic spectra are visible at large dis- 
tances (z > 0.1), because soft spectra more eas- 
ily fall below the current instrumental sensitivity 
limits. Another explanation for the prevalent hard 
spectra at large z is an overcorrection of the EBL 
attenuation effects. None of the nearby sources, for 
which no strong EBL modifications apply, show Y 
much smaller than 2.0. Additionally, the detected 
number of objects with soft spectra increased sub- 
stantially since 2002. 

While there is no obvious correlation between M, 
and the VHE 7-ray luminosity, it might be that 
the current data populate only a certain area in the 
M. — r plane. Owing to the large uncertainties of 
the M. determination and the still poor statistics, 
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Figure 3: Correlations of T, v 1 L 1 with M„ z. PG 1553 at assumed z — 0.1 and z = 0.3. 



the future will have to show if such trends are real. 
Perhaps the VHE 7-ray emission is more sensitive 
to the BH spin, the accretion rate or, more impor- 
tantly, of the acceleration environment rather than 
the BH mass. Also results on timing properties 
(see below) support such claims. 

The Case of PG 1553+113 

PG 1553 is a recently discovered TeV blazar [14] 
with unknown distance. With increasing redshift 
z, the intrinsic luminosity has to increase stronger 
than quadratic (v 1 L 1 > 4ir ■ d 2 F) due to EBL 
absorption as to sustain the measured VHE flux 
(Fig. 4). We assume here that PG 1553+113 is an 
"off the shelf" blazar, i.e. with no extraordinarily 
high L. This assumption is difficult to quantify, 
but when translating it into the limit that L is not 
more than 30 times higher than the highest lumi- 
nosities observed, one obtains z < 0.48 (2a limit). 
An extreme luminosity 1000-times higher yields a 
limit of z < 0.68. Among the extreme BL Lac 
objects we find (V 7 L 7 ) max < 10 45 4 erg sr _1 s _1 
for Mkn 501 in the flare state. These limits do not 
only depend on a good knowledge of the EBL at- 
tenuation over a wide range in redshift, but also 
on an assumed reasonable maximum VHE blazar 
luminosity that is strongly dependent on the jet 
Doppler factor 6. In any case, either a strikingly 
high luminosity or a very high S is needed should 
PG 1553+1 13 be more distant than z > 0.35. Pre- 
sumably such very extreme objects are so rare that 
a sufficiently large volume had to be probed to find 
one of them. 
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Figure 4: Luminosity evolution for PG 1553+113 
at different assumed source distances. 

X-ray Duty Cycle and VHE Variability 
Time Scale 

Following a method described in [21] we deter- 
mine the time fraction at which the 2-10 keV flux 
exceeds 50% of the average flux ("duty cycle", 
DC). In addition we require this deviation to be 
significant (S > 5er). Note the outstanding DC of 
Mkn 421. Supporting the claim that variability is a 
defining property of BL Lacs, a flat distribution of 
the DC in v 1 L 1 is found (Fig. 5). A previous study 
[21] including Mkn 421, Mkn 501, 1ES 2344, 
1H 1426 and 1ES 1959 only had found indica- 
tions for an anticorrelation of DC and M,, which 
in our enlarged sample is weakened mainly by the 
recently discovered sources 1ES 2356, PKS 0548, 
and BL Lac. 

Turning to the minimum VHE variability 
timescales r, these do not scale with M.. This 
implies that flares originate from a much smaller 
region than the BH radius and (more importantly) 
that the BH properties do not influence the emis- 
sion process too much, but the jet environment 
may be more important. Note that, in spite of the 
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expected scaling behavior the TeV blazars hosting 
the more massive BH, Mkn 501 and Mkn 421, 
seem to exhibit the smallest r. This, however, may 
be a selection effect caused (1) by their proximity, 
and (2) by instrumental sensitivity, as small r 
measurements require strong sources. The latter 
also disables strong claims about r-luminosity cor- 
relations yet, and all r values are to be understood 
as upper limits. 
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Figure 5: Correlations of X-ray duty cycle and 
VHE variability scales with M. and u 7 L T 

Conclusions and Outlook 

The observation of VHE blazars has started to 
become less biased: Not only blazars with hard 
spectra or in a flaring state are now detected, 
but a much higher dynamical range of VHE 7 
emission levels and states is probed, flare statis- 
tics studies (e.g. [22]) are within reach, and 
generic blazar properties start to become accessi- 
ble. Thus the era of VHE blazar astronomy has 
been entered — astronomy being understood as the 
study of generic properties of a given class of ob- 
jects. 
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